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1. Introduction ii /k23 

In studying the transmission of heat in fluids subject to 
supercritical pressure, we observe, for relatively low thermal 
flow values, a major Increase in the heat transmission coefficient 
when the average temperature of the fluid mixture in turbulent 
flow within a tube approaches its pseudocrltical temperature. 

There are differences of opinion in the attempt to explain 
the mechanism which occurs in this process of heat transmission. 

Some researchers, such as Knapp and Sabersky [1]', Yamagata et. al. 
[2], Goldmann [3] and Shitsman [4], assume that to explain this 
singular behavior of the heat transmission coefficient in the 
neighborhood of the pseudocrltical temperature one must resort to 
the existence of a phenomenon similar to that of nucleated or 
pelicular boiling in subcritical fluids. This leads to the theory 
of pseudoboiling in supercritical fluids for whose explanation it 
would perhaps be necessary to postulate the transient coexistence 
of two phases at supercritical pressures. 

Other researchers, such as Bourke et al. [5], Tanaka et al. 

[6] and Swenson, Carver and Kakarala [73, assume that this singu- 
larity in the neighborhood of the pseudocrltical temperature can 
be explained as being caused merely by a major variation of the 
thermodynamic and transfer properties at pseudocrltical temper- 
ature . 

Although it may be difficult to select one or another theory, 
taking into account the experimental observations of the phen- 
omenon, this author agrees with the researchers who defend the 

* Numbers in the margin indicate pagination in the foreign text. 
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second theory, l.e., he is of the opinion that the major variation 
of the thermodynamic and transfer properties at pseudocritical 
temperature, conveniently taken into account in the function 
performed by the spins can then explain the singularity of the 
membrane 1 s coefficient in its neighborhood. 

Various studies, both experimental and theoretical, have 
already been published on this subject. For the case of forced 
turbulent flow of a supercritical fluid within a tube, the theo- 
retical analyses by Deissler [8], Wllderecht and Sonnernann, ex- 
tended to the case of supercritical fluid by Hess and Kunz [9], 
Goldmann [3], Tanaka et al. [6] and Shiralkar and Griffith [10] 
are well known. These analyses take into account the shear stress 
and the thermal flow, both Invariable, and with a linear variation 
in the radial direction, which would be subject to discussion for 
the case of supercritical fluid which represents a high degree of 
property variation in the radial direction in nonlsothermal flow. 

All the theoretical analyses mentioned also assume the equality 
of the thermal spin diffusibilities and of the momentum, e M » e t /k25 
and differ essentially from one another only in the ratios employed 
to determine the value of the spin diffuslbility of the momentum. 

In the theoretical analysis of this study the supercritical 
fluid flow is being considered within a smooth cylindrical tube in 
a turbulent regime with the profile of velocity fully developed. 

For the sake of simple presentation, the cases of constant and 
axi symmetrical thermal flow are being considered. 

In order to approach real conditions as much as possible, 
the following is being considered in this study: 

a - a real variation of a thermal flow in the radial direction 
obtained by the application of the first principle of thermodynamics; 

b - a linear variation, in the radial direction, of the shear 
stress; 

c - a ratio between the thermal spin diffusibilities and the 
momentum different from unity (the turbulent Prandtl number dif- 
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fers from unity, a function of the point considered, for each 
set of values of the parameters which influence this process. 


2. Basic Equations 

The equations of momentum and energy. In the forms in- 
dicated below, are the fundamental equations which, integrated 
for each section of the tube, will yield the distributions of 
velocity and temperature: 


t * t> < v + c H 


du 

dy 


o 


p <= 


-( a + c T ) 


dT 

dy 


( 1 ) 

( 2 ) 


where t = shear stress; p = density; v = kinetic viscosity; 

= spin diffusibility of the momentum; a = thermal diffusibil- 
ity; e T = spin thermal diffusibility; y = velocity component in 
the principal flow direction (x); T = temperature; y = coordin- / i t26 
ate perpendicular to the tube surface, with the origin at this 
surface. 


To permit considering the variation of properties that occcur 
in the supercritical fluids in the universal distribution of vel- 
ocity, u = f(y T ), the nondimensional velocities and 
distances are defined according to Goldmann [ 3 ], by: 


and 



du 

* 

u 


e 



dy 


(3) 

(*> 


where u* is the shear or attrition velocity, defined by 


u 


* 



(5) 
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If we introduce these expressions in the equation of 
momentum, wc obtain 

-j_ = tl + Jh_ , -dal (6) 

T 0 v dy + 

which formally is the same as that obtained in the case in which 
the properties are considered constant. Thus, the same formal 
expression of the spin diffusivity of the momentum obtained for 
the case of fluid flow with constant properties can be employed, 
provided expressions (3) and (*l) are being used. 


The spin diffusivity of the momentum used is that developed 
by Spalding [11]: 



1,0 + 0, 04432 [" exp (0, 4u + ) -1, 0-0, 4 


u 4... «Mu+) 2 

2! 


3 ■ 


I'O 


( 7 ) 


The thezvnal spin diffusivity used is that obtained by Yamane /H27 
[12] by means of a modification in the theory of the frandtl 
mixture length and is based on the calculation of transfer effi- 
ciency of momentum and energy: 



6 n M *r c m 1 *r 

1, 0+exp{ [2-0+0, GN ( “ r — ) ]} 

Np ( _2L) n M 

* V* ' >1 * 


[(0,40972 ) Z +0 , 088 0,40972 


(8) 


where 


n M -0,5{ [(0,40972 


M ) 2 +0,8888 2 -0,40972 


v 


3. Results 

Considering a linear distribution of the shear stress and 
a thermal flow distribution given by the first principle of 



thermodynamics, the fundamental equations and the differential 
equations which define the characteristic parameters are inte- 
grated by the numerical method of Runge-Kutta-Gill and furnish 
the velocity distribution, the temperature distribution, the 
average temperature of the mixture (based on the average enthalpy 
of the mixture), and the membrane coefficient for given values 
of thermal flow on the surface (q Q ) and of the discharge (tfi). 

The Runge-Kutta-Gill method is applied in this case by means 
of an iteration process with respect to the shear stress on the 
surface whose value is a discharge function; and due to the var- 
iation of the thermodynamic and transfer properties as a function 
of temperature, the equations are simultaneously integrated. 

In this manner, for the case of carbon dioxide whose pro- 
perties are reasonably well known, we can obtain graphs such as 
those furnished by the velocity distribution, the temperature 
distribution and the membrane coefficient for given values of 
the characteristic parameters. 

Fig„ 1 represents the comparative graphs of the results ob- 
tained from this theoretical analysis with the experimental re- 
sults and of the theoretical analysis performed by Tanaka et al. 

[6] and shows good agreement with the experimental results. Par- /^23 
tleularly Fig. 2, which locates the elements of the fluid at a 
pseudocritical temperature is the one which is of interest to the 
thermal analysis of this study. 


'4. 


Thermal Analysis 
From the equation of energy 


at; 


t -t = / - 
» 


total o c o J k+p e ^ 


dy 


(2), we obtain 
- {at) scl +ut) ca +(at) ct 


(9) 


where y = the value of y at the center of the tube; SCL refers 29 

Q ^ 

to the laminar sublayer (y between 0 and 5); CA refers to the 
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Fig. 1. Comparison between the theoretical 
curves of the heat transfer coefficient and 
the experimental points. 

theoretical curve by author. 

theorec'cal curve obtained by Tanaka e 

experimentation 




buffer layer (y + between 5 and 30); CT refers to the turbulent 
layer (y >30) and k = thermomolecular conductivity. 

The heat transfer coefficient in the case of turbulent 
flow may be calculated approximately as follows: 



where T = temperature of the internal tube surface; T = aver- 
o m 

age temperature of the mixture; T = temperature at the center 

c 

of the tube and h = coefficient of heat transfer. 

Therefore, the membrane coefficient will have a high value 
when (A^)fcotal SIBa ^> l.e. when (k + p CpE^) is large. This 
will happen when the temperature at certain points in the fluid 
is in the zone of the pseudocrltical temperature, because at 
that temperature k and pc p will have maximal values. 

However, so that the effect of this major increase in val- 
ues of the properties at pseudocrltical temperature may become 
noticeable, causing a reduction in the thermal resistance, it 
is necessary for the fluid element which is at the pseudocrltical 
temperature to be located within the buffer zone, i.e. in the in- 

4 . 

terval of the nondimensional distance, y , between 5.0 and 30.0. 

If the fluid point where the temperature is pseudocrltical 
is located in the turbulent layer, the effect of the increase 
of those properties, K and pc , will be small because in this layer, 
the thermal resistance is already much smaller with respect to 
the other layers due to the high value of the spin thermal dif- /43Q 
fusivity e^. Therefore, in this case, the increase in the mem- 
brane coefficient due to the increase in the property values will 
not be significant* 

If the point whose temperature is pseudocrltical is located 
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in the laminar sublayer, the contribution to the increase in the 
membrane coefficient will be only k because (p c„ in this 

“ pi 

layer is negligible duo to the value of the dlffuslvlty e^, And 
since the increase of k is relatively small compared to the spec- 
ific heat, the effect of the variation of properties upon the in- 
crease of the membrane coefficient will also not be significant, 

If on the other hand, the fluid element at pseudocrltical 
temperature is located in the buffer zone where the terms k and 
p 0 ^ are of the same order of magnitude and where (AT) CA con- 
stitutes the major portion of (^^otal’ the simultaneous in- 
crease of k and p c at the pseudocrltical will yield a consider- 
able decrease in (AT)^ 0 ^ In this manner, greater values of 
the membrane coefficient; will be obtained resulting In the forma- 
tion of those characteristic peaks of the coefficient, shown in 
Pigs. 1 and 2. 

Pig. 2 confirms the analysis performed above. We show the 
curve of the points y for which the temperature is the pseudo- 
critical temperature as a function of the average temperature of 
the mixture for the conditions indicated in the figure. It can 
be verified that the peak of the membrane coefficient occurs 
when the point of the fluid which is at the pseudocrltical temper- 
ature is located in the interval 5 < y <30 corresponding to the 
buffer layer. 

5, Conclusions /A31 

The theoretical analysis of the process of heat transfer of 
a supercritical fluid, flowing in a turbulent regime within a 
tube, performed in this study, leads to results which present 
good agreement with the experimental results in the region of 
the peak characteristic of the membrane coefficient. This analysis 
permits an exolana.tion of the formation of the characteristic peak 
and justifies the necessity of the fluid element at pseudocrltical 
temperature to be located in the buffer layer for a major increase 
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in the heat transfer coefficient to occur. By similar reasoning 
to that of the previous section, the thermal analysis of the 
mechanism of the heat transfer process also permits us to explain 
th^ following: 

a - the formation of the peak of the heat transfer coefficient 
at an average mixture less than in the case of heating; 

b - a reduction in the value of the peak of the heat trans- 
fer coefficient as a function of the Increase in the absolute 
value of the thermal flow when the discharge is maintained con- 
stant; 

c - obtaining a higher value for the peak of the coefficient 
in case of cooling than that for heating, maintaining the absolute 
value of the thermal flow and the discharge constant. 
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